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We report on dynamic properties of a simple model microswimmer composed of three spheres 
and propelling itself in a viscous fluid by spinning motion of the spheres under zero net torque 
. . . constraint. At a fixed temperature and increasing the spinning frequency, the swimmer demonstrates 

' a transition from dissipation-dominated to a pumping-dominated motion regime characterized by 

' negative effective friction coefficient. In the limit of high frequencies, the diffusion of the swimmer 

, can be described by a model of an active particle with constant velocity. 

(N ■ 

Attention of physics community to the problem of swimming at low Reynolds numbers, which is relevant for the 
I world of microorganisms, was attracted by Edward Purcell in seventies He formulated the basic principles of 
' self-propulsion and suggested a variety of simple model microswimmers that would propel themselves in the Stokesian 
regime using non-reciprocal cyclic moves. The most famous of them, three-link swimmer, was solved analytically only 
recently Q]- As the Stokesian regime is characterized by the absence of time in the flow equations, the description of 
self-propulsion reduces to a purely geometrical problem of transformation of the microswimmer's body shape. The 
problem was solved for various nearly spherical objects, whose surface is deformed by a wave-like perturbation in the 
manner of ciliated microorganisms [3, US, S 01 ■ number of other simple models performing one- or two-dimensional 
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^ ' non-reciprocal moves as well as their swimming performance were discussed recently 0, [l^ . 
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Recent advances in micromanipulation techniques made it possible to construct artificial swimmers mimicking the 



' bacterial and protozoal self-propulsion mechanisms. Most of these machines, however, are supposed to be driven by 
an external fields rather than ATP hydrolysis. The first working device imitating the flagellum beating and driven by 
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^ ' an external magnetic field was reported recently ll| . Realistic implementations of DNA-based nanomachines using 
the ratchet principle were also suggested [l^ . ITsj . Other directions in development of self-propelling micromachines 
O is related to use of anisotropic environments, active surfaces or chemical reactions, whose mechanical response has 
inherent asymmetry 14, 15, 1^, 17, H, [l^l- The focus of these publications is the propulsion mechanism as such 
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^ ' and the dynamics of a swimmer on long timescales is usually not addressed. One should note, however, that on the 
^ , nanometer and micrometer scale the thermal fluctuations are expected to compete with the propulsion mechanisms 
(SJ ' and, therefore, the interplay between the swimming and dissipation processes is of great interest |2ll . 
On ■ In this work, we report on general dynamic properties of a microswimmer at finite temperatures. 

For our study we chose a simple model swimmer, consisting of three spheres with their centers comprising an 



equilateral triangle. The distances between the spheres are fixed. To impose propulsion, we make the spheres spin via 



On ■ applied constant torque. We arrange the torque directions as shown in Fig. la so that the net torque is always zero. 



This constraint mimics the rotation due to internal degrees of freedom, as it would be in case of a microorganism 
swimming. At the same time, the algorithm allows us to keep control over the amount of supplied energy. The 
dynamics of the three-ball animals was modeled numerically using a hybrid molecular dynamics/ Lattice Boltzmann 



' (L-B) simulation method [27|, |28|. The spheres were modeled with the raspberry setup (a tethered network of point 
^ 1; particles wrapped around a big Lennard-Jones (LJ) bead) [11] with the radius i? = 3, while the links by the massless 
^ harmonic bonds of length L = 10 of high spring constant between the ball centers. The length units are defined by 
the LB cell size. Most simulations were performed in a cubic box of side 40 units with periodic boundary conditions. 
The dynamic viscosity of the fluid in the LJ units was set to 77 = 2.55 if not specified otherwise, and the temperature 
to r = 5. To simulate the action of finite temperature, a fluctuating stress tensor technique was used and in addition 
stochastic random force was applied to each of the beads composing the swimmer j^^]. The velocity of each bead v 
is then calculated from the following Langevin equation 

m^ = -C(v-u(r))+f(i) + Fext+Fint (1) 
at 

where m is the bead mass, C, the friction coefficient, u(r)) the interpolated fiuid velocity at the particle position r, f{t) is 
a Gaussian white noise force with zero mean, whose strength is given via the standard fiuctuation-dissipation theorem 
in the absence of the external force, Fint the interaction force between different beads comprising the raspberry (see 
Ref. [2^ for details), and Fgxt the "external" force driving the sphere rotation. The magnitude of the external force 
was chosen to produce the desired spinning velocity U3— T/ {Sr/TrR^), where the net torque T = J^iPi ■ F^ is a sum of 





FIG. 1: The construction scheme of the three-baU swimmer: The baUs are spinning around the in-plane axes so that the angular 
velocities cj2+'i'3 = 0. A simulation snapshot on the right hand side illustrates the raspberry model implementation of 
the swimmer. 

the moments acting on each surface bead of the sphere, p is the distance from the rotation axis to the ith bead. The 
energy influx is then given by a sum e = F • over all beads constituting the spheres or simply e = T-a;. 

Motion of the swimmer at zero temperature can be easily understood in the limit of small spheres and long links. 
Each sphere can be then represented by a rotlet [2§| 
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v(r) — —oj X r 



(2) 



We have to find the flow field produced by two small spheres at the location of the third one and take just sum of 
their the components normal to the plane containing the three sphere centers, as the other components cancel. To 
have a zero net force acting on each of the spheres, it requires the sphere to move normally to the triangle plane with 
the velocity 
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We tested this prediction against simulation results for a swimmer with R — 3 and L — 10 (see Fig. [5]). Despite 
the large ball size as compared to the link length, the equation gives a very good velocity estimate with relative 
deviation from the simulation data less than 4% up to w = 2. The result given by Eq. ([3]) can be also be compared 
to the propulsion speed of a twirling ring, which corresponds to a generalization of the three-ball animal to the case 
of inflnite number of balls filling a full circle of radius r 12, isj : 
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The propulsion efficiency of this swimmer can be estimated as a relation of the energy dissipated in the forward 
thrust to the energy dissipated in the sphere spinning, which gives 



/3 = 



Fftvo 
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(5) 



The steady state thrust force Fft is equal to the total friction force Ff due to forward motion of the three-ball 
swimmer, which is related to the corresponding mobility as Ff — —vo/fis. If the hydrodynamic interactions between 
the spheres are taken into account in the stokeslet approximation (the mobility of a pair of stokeslets to the first order 
in r^^ is = ^(1 + jy)), we get for the triplet 
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where /zi = {GirriR) ^ for the motion transversal to the triangle plane. The torque of the friction force exerted by 
a spinning sphere of radius R is given by T = STrrjR^. If we neglect hydrodynamic interactions (due to fast decay 
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FIG. 2: Self-propulsion velocity of the three-ball swimmer with R = 3 and L = 10 as a function of the ball spinning frequency. 
The solid line was obtained using Eq. [S] 

of the the flow field magnitude for a rotlet, (R/r)^^, it remains a good approximation for (R/L) < 1), the energy 
dissipated just due to independent spinning of three spheres becomes Stt/jR^lo. The efficiency is then (3 ~ 0.5(i?/L)*, 
which reaches its maximum of about 6% for contacting spheres. 

We now study the motion of the swimmer at finite temperatures. Typical two-dimensional projections of the 
swimmer's trajectories observed at different ball spinning frequencies are plotted in Fig. [3) One can see that the 
character of the motion is changing significantly on increasing the frequency: the trajectories become more persistent 
with more extended straight segments and less abrupt changes of the direction. As a result, the swimmer with a 
higher lu travels over longer distances within the same time as compared with a less active one. 

The mean-square displacements (MSD) of the swimmer's center of mass are shown in Fig. for four different 
spinning frequencies and three different temperatures. We see immediately, that the higher spinning frequencies 
correspond to a higher MSD, i.e. to the higher diffusivity. The shape of the MSD curves clearly indicates the 
existence of two distinct dynamic regimes: the driven motion at t < 500 where (Ar^) cx t'^ and the diffusive motion at 
t > 1000 where (Ar^) cx t. In the case of high propulsion speeds, the transition from the deterministic driven motion to 
a random walk happens due to rotational diffusion of the swimmer. The "ballistic part" of the MSD coincides for the 
swimmers with the same frequency but different temperatures. In contrast, the random walk part is thermosensitive: 
higher temperature leads to a lower diffusivity. The transition between the two regimes happens on the characteristic 
timescale of rotational diffusion. The diffusion coefficient measured from the asymptotic behaviour of the swimmer 
at long times, I?cff, is plotted in Fig. [4Jd. We see that the measured diffusion coefficient grows nonlinearly with the 
drift velocity (or the spinning frequency as vq cx w). 

A simplest model describing the diffusion of an active swimmer moving with a constant velocity vq, whose direction of 
motion changes due to rotational diffusion with a diffusion constant Dr, gives the diffusion coefficient D^ff = 2vq/ Dr- 
This expression follows immediately from a reformulation of the known result for the end-to-end distance of a persistent 
polymer chain, which is (Ar^) = 2lpL, Ip being the persistence length and L the contour length of the chain. In the 
swimmer variables, the persistence length of the trajectory is Ip = tvq, where r is the rotational correlation time, 
and its contour length is simply proportional to time: L = tvQ. Then, D^ft = (Ar^)/6i = Vqt/3 [20]. The time 
r is defined via the orientational correlation function of the swimmer (n(0) • n(t)), where n is the normal vector to 
the plane of the triangle connecting the sphere centers. In case of random uncorrelated fluctuations, the correlation 
function decays at short times exponentially with time as (n(0) • n{t)) = e"'/"^. In case of an anisotropic rotation, one 
can expect that the decay is governed by the smallest correlation time, i.e. by the characteristic time of rotation in 
the direction with the smallest friction torque. The shortest orientational relaxation time expected for this swimmer 
corresponds to a rotation around the in-plane axis passing through the center of one ball and through the middle of 
the spring connecting the opposite ones. With a first order correction for hydrodynamic interaction between the balls 
we get 
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FIG. 3: Two-dimensional XY projection of the trajectories of the swimmer's center of mass position obtained with different 
ball spinning frequencies at T = 5. The insert shows the magnified view of the initial part of the trajectories. 

This formula gives r — 2250, 1350, and 675 for T ~ 3, 5, and 10, respectively. This values are close to what we 
observe in simulation at low to. If we use this correlation time to predict the translational self-diffusion coefficient, the 
agreement with the measured data is convincing (the solid curves in Fig. [Hd). We also tested our assumption regarding 
the orientational correlation time by direct measurement of the correlation function of the swimmer trajectories. From 
these measurements at low spinning frequencies we get r = 2300 ± 200, 1300 ± 200 and 700 ± 100 for T = 3, 5, 
and 10, respectively. These values are in good agreement with the estimated correlation times. In what concerns 
the assumption of constant velocity along the trajectory, we can check the swimmer's velocity distribution. The 
characteristic distributions of the velocity component perpendicular to the triangle plane (i.e. parallel to the thrust 
force) for T = 3 are shown in FigO While at w = 0.24 the variance is comparable to the mean velocity, at the 
highest frequency, w = 1.2, it does not exceed 25% of the mean. So, we see that the swimmer's velocity is narrowly 
distributed around vq and that also this assumption is justified. Thus, the swimmer's motion, especially at the high 
ball spinning frequencies (or energy influx rates) fits well the model of an active Brownian particle [ill . [23| , whose 
velocity is fairly constant in the magnitude while the direction is subject to thermal fluctuations. In other properties 
we observe more analogies to the active Brownian particle model. For example, the existence of a finite preferential 
velocity in our model is also reflected in the velocity distribution function shown in Fig. [5) Instead of a maximum at 
V = characteristic for systems in thermal equilibrium (the Maxwell distribution) it possesses two maxima at ±Vp 
and a minimum at u = 0. 

At Vq > (u^(w = 0)) we expect that the character of the swimmer's motion changes. As it was mentioned above, 
the characteristic velocity vq is a function of the energy influx rate, while the variance is fixed by the temperature and 
the swimmer mass: (t;^) = iksT/M , and characterizes the energy dissipation rate. Thus, the transition is expected at 
Lj = kBTL^/{R^M). The threshold values of the frequency are w « 0.3, 0.5 and 1.0 (or = 0.13, 0.23, and 0.45) for 
T = 3, 5 and 10, respectively. At the higher frequencies, other signatures of the driven regime can be observed, such 
as negative friction coefhcient [23| . The friction coefficient can be extracted from the relation between the swimmer 
acceleration a and its instant velocity v . We calculated these values by averaging over three velocity/acceleration 
components. The result for different spinning frequencies is shown in Fig. [6l top. As there is no other external 
forces acting of the swimmer, the acceleration is resulting from the viscous friction a — Ff^/M = —j{v)v/M. Here, 
the effective friction coefficient "f{v) is a function of the instant velocity. The a{v) curves for the lowest spinning 
frequency are fairly linear at low velocities, which corresponds to a constant positive friction coefficient. The curves 
bend toward smaller friction at the higher frequencies, and at the highest one all the curves cross zero and display 
a region of positive acceleration. All the a(v) curves show two linear regions: one at u — > and another one at the 
higher velocities. The corresponding nonlinear friction coefficient is plotted in Fig. [S] bottom. There we indeed notice 
saturation both at low velocities and at high v. While the limiting values 7(1; — > 0) shows strong dependence on 
the energy influx rate, the high velocity values seem to tend to a temperature independent limit, which only slightly 
decreases with the ball spinning frequency. We also notice that the parts of the curves with the acceleration parallel 
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FIG. 4: Top: Mean-square displacement of the swimmer as a function of the ball spinning frequency and temperature. Bottom: 
Long-time diffusion coefficient of the swimmer at different temperatures measured from the long-time asymptotic behavior of 
the mean-square displacement. The solid lines are calculated as DcB = 2vo/Dr (see the text for further details). 




FIG. 5: Distribution function of the swimmer's velocity component in the direction of the thrust force at different spinning 
frequencies and temperature T = 3. 
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FIG. 6: Top: Swimmer acceleration as a function of its instant velocity at different ball spinning frequencies (marked on the 
curves) and temperatures. Bottom: The corresponding effective friction coefficient as a function of the swimmer's velocity. 



to the instant velocity correspond to 7(w) < 0. One can see that the transition into the driven regime with negative 
7 occurs on increasing lo. At T = 3 and T = 5 the negative friction appears at w = 0.9 and 1.2, at T = 10 only 
at w = 1.2, in agreement with the above estimates. The values Vp, at which the friction coefficient turns zero reflect 
therefore the preferential stationary velocity, at which the viscous resistance of the medium is exactly balanced by 
the swimmer's forward thrust. 

We would like to note also that the behaviour of the friction coefficient qualitatively resembles the velocity depen- 
dence predicted within the active Brownian particle model . In the adiabatic approximation the effective friction 
for a particle with an internal energy depot should vary as 7(w) = 7o ~ fR^' ^^^^''^ the coefficient a is proportional 
to the energy influx rate. If applied to our case, the 70 would be the Stokesian viscous friction experienced by the 
swimmer. The negative values of the nonlinear friction coefficient can be observed provided that 7(1;) turns negative 
at finite velocities, which requires that a > 70(6 + v^). We can see from FiglHlat 0.24 < w < 0.9 that the shape of the 
friction curves 7(w) shows a similar behaviour. 

Finally, we would like to note that the parameters of this system were chosen to demonstrate the wide range 
of dynamic properties. In case of a micrometer size living swimmers, the role of thermal fluctuations might be 
signiflcantly smaller. Living microswimmers moving at the velocities of lO/xm/s would rather demonstrate the driven 
motion regime. Still, in presence of a stochasticity caused by internal mechanisms of the animal or spatial/ temporal 
inhomogencities of the medium, its competition with the deterministic drive would lead to qualitatively similar results. 
On the other hand, slow and small swimmers like the artificial ones, suggested in literature, would demonstrate many 
features described in this work even with thermal fluctuations. Another parameter that differentiates the dynamics 
studied in the present work from the biological microorganisms is the relatively high Reynolds number that reaches 
Rsmaa; ~ 3 for the highest velocities considered. So, the role of inertial effects is considerably overestimated as 
compared to the microworld. Our tests, however, showed that none of the reported properties can be attributed solely 
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to inertia and we expect our predictions to be valid also for the typical microscopic objects. 

We studied dynamics properties of a self-propelling microswimmer engine subjected to thermal fluctuations. Our 
model for the first time treats explicitly both the propulsion mechanics and fluctuations and demonstrates a rich 
interplay of these factors. Qualitatively, many of the observed properties (such as the non-linear velocity-dependent 
friction coefficient, characteristic crater-like velocity distributions at high energy influx rates, etc.) resemble those of 
the model of active Brownian particle [ll], Ull- The differences related to the fluctuations of the propulsive effort, the 
broad velocity distributions, three-dimensional motion, and anisotropic friction, result in general in a more complex 
dynamic behaviour, which we hope will stimulate further theoretical effort in this field. 

We are grateful R. Thaokar, W. Ebeling, L. Schimansky-Geier, U. Erdmann, B. U. Felderhof, J.-F. Joanny for 
stimulating discussions. 
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